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Abstract 

We describe the CGMY and Meixner processes as time changed Brow- 
nian motions. The CGMY uses a time change absolutely continuous with 
respect to the one-sided stable (Y/2) subordinator while the Meixner time 
change is absolutely continuous with respect to the one sided stable (1/2) 
subordinator. The required time changes may be generated by simulating 
the requisite one-sided stable subordinator and throwing away some of 
the jumps as described in Rosinski (2001). 



1 Introduction 

Levy processes are increasingly being used to model the local motion of asset 
returns, permitting the use of distributions that are both skewed and capable of 
matching the high levels of kurtosis observed in factors driving equity returns. 
By way of examples we cite the normal inverse Gaussian process (Barndorff- 
Nielsen (1998)), the hyperbolic process (Eberlein, Keller and Prause (1998)), 
and the variance gamma process (Madan, Carr and Chang (1998)). For the 
valuation of structured equity products the importance of skewness is well rec- 
ognized and has led to the development of local Levy processes (See Carr, Ge- 
man, Madan and Yor (2004)) that preserve skews in forward implied volatility 
curves. It is also understood from the steepness of implied volatility curves that 



1 



tail events have significantly higher prices than those implied by a Gaussian 
distribution with the consequence that pricing distributions display high levels 
of excess kurtosis. 

On a single asset one may simulate the Levy process calibrated to the prices 
of vanilla options to value equity structured products written on a single un- 
derlier. Such a simulation (See Rosinski (2001)) may approximate the small 
jumps using a diffusion process with the large jumps simulated as a compound 
Poisson process where one uses the normalized large jump Levy measure as 
the density of jump magnitudes with the integral of the Levy measure over the 
large jumps serving as the jump arrival rate. However, increasingly one sees 
multiasset structures being traded and this requires a modeling of asset correla- 
tions. Given marginal Levy processes one could accomodate correlations if one 
can represent the Levy process as time changed Brownian motion. In this case 
wc correlate the simulated processes by correlating the Brownian motions while 
preserving the independent time changes for each of the marginal underliers. 

It is therefore useful to have representations of Levy processes as time 
changed Brownian motions. For some Levy processes, like the variance gamma 
process or the normal inverse Gaussian process, these are known by construction 
of the Levy process via such a representation. For other Levy processes, like 
the CGMY process (Carr, Geman, Madan and Yor (2002), see also Koponen 
(1995), Boyarchenko and Levendorskii (1999, 2000)) or the Meixner process 
(Schoutens and Teugcls (1998) see also Gregelionis (1999), Schoutens (2000), 
and Pitman and Yor (2003)), the process is defined directly by its Levy mea- 
sure and it is not clear a priori whether the processes can be represented as 
time changed Brownian motions. With a view to enhancing the applicability 
of these processes, particularly with respect to multiasset structured products, 
we develop the representations of these processes as time changed Brownian 
motions. 

Section 2 presents for completeness, some preliminary results on Levy pro- 
cesses that we employ in the subsequent development. In section 3 we develop 
the CGMY process as a time changed Brownian motion with drift, where the 
law of the time change is absolutely continuous over finite time intervals with 
respect to that of the one sided stable Y/2 subordinator. The simulation of 
CGMY as time changed Brownian motion is described in section 3. Section 4 
develops the time change for the Meixner process as absolutely continuous with 
respect to the one-sided stable 1/2 subordinator. Simulation strategies for the 
Meixner process based on these representations are described in Section 5. Sec- 
tion 6 reports on the simulation results using chi-squared goodness of fit tests. 
Section 7 concludes. 

2 Preliminary results on Levy processes 

We present three results from the theory of Levy processes that we make critical 
use of in our subsequent development. The first result relates the Levy measure 
of a process obtained on subordinating a Brownian motion to the Levy measure 
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of the subordinator. The second result establishes a criterion for the absolute 
continuity of a subordinator with respect to another subordinator. The third 
result presents the detailed relationship between the standard presentation of 
the characteristic function of a two sided jump and one-sided jump stable Levy 
process and its Levy measure. These are presented in three short subsections. 

2.1 Levy measure of a subordinated Brownian motion 

Suppose the Levy process X(i) is obtained by subordinating Brownian motion 
with drift (i.e. the process 9u + W(u), for (W(u),u > 0) a Brownian motion) 
by an independent subordinator Y(t) with Levy measure v{dy). Then applying 
Sato (1999) theorem 30.1 we get that the Levy measure of the process X(t) is 
given by fi(dx) where 

r°° i (*-o«) 2 

mdx) — dx I v(dy)—==e 2 » . (1) 
Jo V 2 7n/ 

2.2 Absolute Continuity Criterion for subordinators 

Suppose we have two subordinators T A = {T A {t),t > 0),T B = (T B {t),t > 0). 
The law of the subordinator Ta is absolutely continuous with respect to the 
subordinator Tb, on finite time intervals, just if there exists a function f(t) such 
that the Levy measures VA(dt), VB{dt) for the processes Ta and Tb respectively 
are related by 

v A (dt) = f(t)u B (dt) (2) 
and furthermore, (Sato (1999) Theorem 33.1) 

J v B {dt) (yjtfj- if <oo. (3) 

2.3 Stable Processes 

The Stable Levy process S(a,a,(3) = (X(t),t > 0) with parameters (<r,a, /3) 
( For details see DuMouchel (1973, 1975), Bertoin (1996), Samorodnitsky and 
Taqqu (1998) Nolan (2001), Ito (2004) ) has a characteristic function in standard 
form 

E[e luX{t) ] = exp(-t*(u)) 
where the characteristic exponent ^(u) is given by 

#(u) = (J a \u\ a (l -i/3sign{u) tan (™)) , a ^ 1 (4) 

= a\u\ ^1 + i(3sign(u)— log(|u|)^ , a = 1. 

The parameters satisfy the restrictions, a > 0,0<a<2 and — 1 < (3 < 1. The 
one sided jump stable processes result when (3=1 and there are only positive 
jumps or [3 = — 1 in which case there are only negative jumps. 
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The Levy density of the stable process is of the form 
and we have that 



a- 



Cp T" C n 



(5) 



(6) 



It remains to express a in terms of the parameters of the Levy measure. In 
the one sided case with only positive jumps we have 



a = 



<*r(f)r(i-f) 



2r(l + a) 

and more generally for the two sided jump case we have 

c P + c n r (f ) r (1 - f ) 

2 r(l + a) 

Conversely, c p and c„ may be computed in terms of /3 and u. 



(7) 



(8) 



3 CGMY as time changed Brownian motion 

We wish to write the CGMY process in the form 

X CGMY (t) = 9Y(t) + W(Y(t)) 

for an increasing time change process given by a subordinator (Y(t),i > 0) 
independent of the Brownian motion (W(s), s > 0) . 
The characteristic function of the CGMY process is 



E[e^{iuX CGM Y{t))] = {4>CGMY{u)f = exp tCT(-Y) 



(M-iu) Y -M Y + 
(G + iuf - G Y 



The complex exponentiation is defined via the complex logarithm with a branch 
cut on the negative real axis with polar coordinate arguments for the complex 
logarithm restricted to the interval ] — ir, +ir]. The CGMY process is defined 
as a pure jump Levy process by its Levy measure 



k 



CGMY 



(x)=C 



cxp(-G|x|). 



ll+Y 



<-x<0 



exp (-Mi) 



rl+Y 



L x>0 



On the other hand we have, in all generality, by conditioning on the time 
change that 



E 



iu{6Y(t)+W{Y(i)) 



= E 
= E 



exp ( iu$Y(t) - ^P-u 2 



exp 



iuO Y(t) 
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Take u(X) to be any solution of 



and in this case 



G + M . / /G-M N 2 



M-iu = — \-i\j2X- 



G + M . I { G - 2 



G + iu = i\2X 



2 y v 2 

It follows that the Laplace transform of the subordinator is 

E[e' XY{t) ] = exp (tCT(-Y) [2r Y cos(r]Y) - M Y - G Y ] ) 
r = V2X + GM 

r/ = arctan 



2A 



>G-M\ 



In the special case of G = M we have 



E[e- XY{t) ] = exp (2iCT(-F) 



Then we have the Laplace transform of the time change subordinator as 
E[e' XY{t) ] = exp (tCT(-Y) [(M - iu(X)) Y — M Y + (G + tu(X)) Y - G Y 
The solutions for u are: 

u = i9± V2X - 2 

where we suppose that 9 2 < 2X. 

We shall see that a good choice for 9 , for sufficiently large A, is 

„ G-M 



(2 A + M 2 ) Y/2 cos (y arctan ( ^ j 1 - 



3.1 The explicit time change for CGMY 

We shall show that the time change subordinator Y(t) associated with the 
CGMY process is absolutely continuous with respect to the one-sided stable 
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Y/2 subordinator and in particular that its Levy measure v{dy) takes the form 

K 



v{dy) 

m 

B 
K 



y 



f(y)dy 



(Bj-Af)y 



e 2 E 
G + M 



g 2 y Ty/2 " 

2 Tl/2 



(9) 



cr 



2r(i + f) 



where 711,71 are two independent gamma variates with unit scale parameters 
and shape parameters Y/2, 1/2 respectively. Further we explicitly evaluate the 
expectation in equation @ in terms of the Hermite functions as follows. 



E 



' B 2 y Ty/2 

e 2 71 / 2 



r(y)r(|) 



where 



I(v,a,X) 



^v—l — ax— An 

x e 



'da; = (2X)-" /2 T(v)h- v 



and h- v (z) is the Hermite function with parameter — v (see e.g Lebedev (1972), 
p 290-291). 



3.2 Determining the time change for CGMY 

For an explicit evaluation of the time change we begin by writing the CGMY 
Levy density in the form 



kcGMY{x) = C- 



Ax — B\x\ 



rl + Y 



-, where: A = 



G-M 



B = 



G + M 



Henceforth, when we encounter a Levy measure n(dx) that is absolutely 
continuous with respect to Lebesgue measure we shall denote its density by 
fi(x). We now employ the result QJ and seek to find a Levy measure of a 
subordinator satisfying 



G- 



Ax — B\x\ 



U+V 



1 (g-"«) 2 

v\dy) „ — e 



I2w\j 



2 ^+9x 



v(dy)— ? ==e 2 v 2 



I2wy 



We set 9 = A and observe that the right choice for 8 is (G — M)/2 as 
remarked earlier, and identify v{dy) such that 



C- 



,-B\x\ 



1 -zL-^JL 

v(dy)—^=e 2 « 2 



rhfy 



(10) 



G 



We now recognize that the Levy measure for the CGMY is (taking C — 

f,/ Y \ 2 i now), that of the symmetric stable Y Levy process with Levy 
T ( 1 +~j 
measure tilted as 

kcGMvix) = e Ax - B ^k staU<Y) {x). 

We also know that 

Xstable(Y)\t) = B Y o( t ) 

where Y°(t) is the one sided stable Y/2 subordinator, independent of the Brow- 
nian motion (B u ) . 
We now write 

X CGM Y(t)=6Y^(t) + W YW(t) 
and we seek to relate the Levy measures and of the processes Y^ and 

From the result we may write 

x 2 

,(0)f- • e 2y 



Mo (a;) = vW(dy) 



r2?ry 

(=c-fly) 2 



2ny 



Hence we must have that 



s/y Jo Vv 

Taking 6 — A, we get: 



Vv Jo Vy 

We now use the well known fact that 
to write 

— % A „ v roc /_> /*oo — ^r( — +^ 



Jo v / 2ttu 3 Jo ' 
By uniqueness of Laplace transforms we get that for every function / 



1+ 



^Hdyf-—f( 1 -)=[ du^£=e-£ [ v(o ){dy) J_ f fl +u 

Vv \vJ Jq \2iru Jo Vv \y ; 
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or equivalcntly that, for every function g : 



Jo VV Jo \/W Jo Ky '^y y \l + uy 

Jo V2iru 3 Jo \l-usj ^1 

Jo u Viw e 2 u Jo (i-su) 2 rz 

V 1_ti 



Hence it is the case that 



k duBe-^^(j^) 
o y/2ir(u(l — uy)Y 
r 1 dvBe-^is^ij^) 
^ Jo ^2tt(v(1 - vfp~ 

In particular we have 

- 1 dvBe'^^"^) v^\^) 



^ (1) (y) = Vy / 

Jo 



^2k{v{1-v)Y 



We now introduce the explicit form of vo(y) f° r °ur case where it is the Levy 
density of the one-sided stable Y/2 subordinator, 

My) 



This gives the representation 

K f 1 dvBe'^^" A2 ) {l-v)^ +1 ) 



Mv) = ~y+T [ 
y 2 Jo 



y/2n(v(l - «)) 



3 



Y+l 

y 2 J i 



K f°° dw „ n^,,,- a 2 \ f w — \ 



y 2 



r dw Bc -%(b* w -a>) 

Ji \J2-kw V w 



KBe-H B2 - A2 ) r dh ^ 

e 2 



y Y 2 X Jo V2ir (1 + h)~ 
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3.2.1 Absolute Continuity relations 

This subsection investigates the absolute continuity relation in general between 
two subordinated processes and the absolute continuity of the subordinators 
as processes. It is easy to show that the laws of the CGMY process and the 
symmetric stable Y process are locally equivalent, i.e. for each t, their laws, as 
restricted to their past a — fields Tt up to time t, are equivalent (from now 
on, as a slight abuse of language, we shall say of 2 such processes, that they 
are equivalent). Now that we have identified these processes as subordinated 
processes, we look for the equivalence in law of the subordinators. Indeed we first 
observe that if the subordinators are equivalent then the subordinated processes 
will be equivalent but the converse may not be true. 
Indeed, consider two subordinators 

T A (t), T fl (t) 

such that the relation J2J) between their Levy measures holds for some function 
f(t) for t > 0. 

We suppose the absolute continuity of Ta with respect to T B or the condition 

©■ 

We also define the subordinated processes 

X B {t) = 0T B (t) 

where (@ u ) is a Brownian motion assumed to be independent of either Ta or 

We have from the result that at the level of Levy measures Ha^b f° r 
Xa, X b 

T 2 

e 2t 



li A {x) = / v A (dt) 
Jo 

POO 

fi B (x) = / u B (dt) 



x 2 



JO V^TTt 

The following then holds as a consequence of ©, for every functional F > 0: 

E [F (T A (s), 8 <t)]=E [F(T B (s), s < t)cf> {T B (s),s < t)] 

where 

*PiW,.so-(^) i 

As a consequence we deduce that, for every G > : 

E[G{XA{s),s<t))=E[G{X B {s), S <i)cj>{T B {s),s<t)] 

Consequently we may write 

E [G (X A (s),s <t)]=E [G (X B (s),s < t) ^(X B (s),s < t)] 
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where 

<p(X B (s), 8 <t)=E [4>(T B (s), s < t)\(X B (s),s 
This implies that we should have 

fi A (dx) = g{x)fi B (dx) 

with 



We want to show that © implies (JUJ. 
Now we have explicitly that 



g(x) - 



Let 



and note that 



We then have 



and 



UD (dt\£- 



Mdt)^- 



g(x) = [ ^ x} (dt)f(t) 



y/rfx)-l=(j \ {x) (dt)f(t)y -1 



(VsM- i) 2 M B (^) = / (Votf)- i) 2 / Mdt) 
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Observe that 



l [x) (dt)f{t) 



v B (dt) 



e 2t 



x 2 

J v B {dt)f{t)^= + j VB {dt) 



v B (dt)- 



'2-Kt 



v B {dt)f{t) 



C 2t 



e 2t 

u B {dt)^= | 



We wish to show that the integral over x of the right hand side is smaller 
than 



v B {dt)f{t) + / v B {dt) - 2 I i/ B (dt)V7(t) 
and this follows provided 



> / *n(<B)v7(t) 



For this consider 



< 



j ' v B {dt)y/Jify = J v B {dt)^W) J dx 

/*/*.<*h/7w(^!)"(5S) 



e 2t 
v/27rf 



2irt 



where we have used Cauchy-Schwarz for fixed x. 
Hence we have 



The result does not go in the other direction as we may take 

v A (dt) = e a {dt) 
v B {dt) = £b{di) 
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for a 7^ b. These are not equivalent subordinators but in this case 



x 2 



H B {x) 




two Levy densities, which in fact are probability densities, so that the cor- 
responding Levy processes which are indeed Compound Poisson, are (locally) 
equivalent. 



3.2.2 Absolute Continuity of the subordinators for CGMY and Sta- 
ble Y/2 

We now establish precisely the absolute continuity relationship between the 
subordinator associated with the CGMY process, and the one sided stable Y/2 
subordinator. 
We note that 

v C GMY{dy) = f(y)v (dy) 



ti ^ -y( B 2 -A 2 ) i-\ Z 00 dh -^-h h 2 

f(y) = e *\ a A )(B^)J^ —e * h 



Y — l 



(1 + hY 



We first check that as B — > for A = we get the expected result that 

For this we let z = B^fy and make the change of variable 

k = z 2 h 

to get 

dk u 



f(y) = e 2 / -=^e 2 



/ 

Jo 



V^z (1 + ^)2 

I 



Y-l 

dk _k (—) 2 
= e 2 I , e 2 . 



I^nk 
dx x 2 

1 



00 dk _k 

e 2 , as z —* 



■■f 

Jo 



For the equivalence of the two subordinators we must check that 



Jo y 2 + v 



1 1 < oo. 
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We break up this quantity in 2 parts dealing with the integral near and oo 
separately. First consider the integral over [l,oo). Here we write 



/ Mdy)f(y)=f -%e-^~ A 2 HBVy)f 

Ji Ji y 2 ^ Jo 

(BVV) I 
Jo 



dk B 2 y, h~ 

-.e ° 



2tt 



(l + hy 



and check that 



dh _B 2 y, 

^=e 2 - 



h 2 



*k {l + h)~ 

is bounded in y. 

Write again B^fy = z, make the change of variable k = z 2 h and observe 
that 



(By/y) 



dh _b^l h 2 
-e 2 



2tt 



(1 + hy 



e 2 — -^- L 



o V2nk 



(1 + ^) 2 



dk _k 
< / e 2 < oo 

/o V27rfc 



We next consider the required integral near 0, or over the interval [0, 1]. We 
have an expression of the form 



B 2 -A 2 



>o y 

_, ,- dh _§\ h I,— 



lo V2^ {l + h) 2 

We now isolate the exponential by writing 



< 



Jo y 2 + v 7 Jo p + v ' 



The exponential term is of order y near zero and hence this first integral is finite. 
For the second one we write 



1 dy 



o y- 



f 1 dy 
Jo y~ +1 

Jo y 2 ^ 
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For the finiteness of this integral we analyse the behavior of (1 — I(y)) near 
zero. For this we analyse I(y) = J(yB 2 ) where 



J(y) = y/y I" 

Jo 

-I 



dh _yh h~^~ 
-.e 2 



2tt 



(l + hy 



dk _k I k 
e 2 $ _ 

/2nk \y 
dx fx 2 

g 2 <J) 

o V^tt V y 

dx x 2 ( x 2 
e 2 $ 



where 



*(0 



i 



Lemma 1 T7ie function is i/ie distribution function of a random variable 

V that can also be realized as the ratio of two independent gamma variates, 
specifically 

where j a is the gamma variate of parameter a. In particular V has finite mo- 
ments of all orders m < 1 and 

Proof. We note that $ is the distribution function of a random variable V 
where for a uniform variate U we have 

/ Y \ 



P(V<0 = P^17< 

= p(u^ < 



1+f 

= p((i + o^<^) 
= p (t/* < c (i - u$ 



so that V is the random variable 



V 



1 - {/-F 
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From the Beta-Gamma algebra we deduce that V is 

V = 



7i 

Consequently V has finite moments for all powers below unity. In particular 
for to < 1 

g[v ro ] = r( / ( y ) m) r(i-m) 

■ 

As a consequence for to = ^ we have that 

r (^±11 



r(f) 

Furthermore we have that as 

i-J(y) = p(\G\<Vvy) 



So the order of convergence of 1 — I(y) = 1 — J(yB 2 ) is always a — ^ and so 

Y 

- <2a=l 

for all Y < 2. The desired absolute continuity result is established. 
We also observe that 

I(y) = J(yB 2 ) = P (|G| > By/vj) = P (J^ > < 1 
3.2.3 A Further analysis of I(y) 

We now write the Levy measure of the CGMY subordinator in the form 

K 



y 



i+- 



-E[e 



-yZ] 



for some random variable Z. 

For a fixed constant B the Levy measure of our subordinator in the sym- 
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metric case is 



R 



KBe 



y+j / 

y 2 j 



dh vB 2 u h 2 
i e 



2tt 



(1 + /^ 



KBe- B i f°° dh y* 2 *^ r 



y 

KBe 



o 



'2tt/i 



Ke~ B ^ 



dk 



Ke 



Jo 



/2ir 
dz 



yBk , . 

e — p(y < k 2 ) 



2ir 



e"P V < 



B 2 y 



G 2 



r + l 



> 



We also know that 



\VB 2 



7i 



with two independent gamma variables. Thus we may write 



R 



Ke 



-B 2 l 



v 2 + 

Ke- B2v i 



-P\li> 



y 2 

K 



+i 



-E 



E 



yB 2 -fY 
G 2 

yB 2 -/i 



exp -- 



G 2 



-+i 



exp -yB 



G 2 



But 



so that we get 



-G 2 = 



R = 



Ke- 



if 



r + l 



-E 



exp -y- 



B 2 ix. 
2"tT 



We now have identified the two Levy measures as 

Kdy 



vo{dy) 



-+i 



and 



v\{dy) 
Z 



(dy)e-^E [exp {-yZ)\ 

B 2 lY/2 

2 7i/ 2 ' 
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3.2.4 Evaluating explicitly the LT of Z 

There is an additional randomness in the simulation if the expectation 

E[e v 2 " 1 /2 ] 

is evaluated by simulation. It is helpful to explicitly evaluate this function, 
begin with 



0..6(A) = E 



exp —A 



76 



Now we have that 



= E 



exp - 



B(b 



0{b,a) 
i f 1 

— / (l-x) a - 1 x b - 1 e--dx 
b , a) Jo 



dy (I 



B(b,a) J 1 y 2 \y 



6-1 



1 

1 - - 

y 



o-l 



-Ay 



dy 

B(b 7 a)J 1 ~y~^ 



(y-ir-'e-^ 



Hence we have that 



~ Xz dz 



B(a,b) Jo (l + z) a + b 
We are interested in the case a = -j, b = | and so we write 



*™ w "*(fi) 



1 / y _ -i X — I 

Y 1 . / du 2 (1 + x) 2 2 e~ 
-, 5) Jo 



Ax 



From Gradshetyn and Ryzhik (1995) (3.38) (7) Page 319 we have 



POO y Y l 

Jo 



Ace 



2-r I ^ ) , 



= 2^r 



(v^A) 



where /^(z) is the Hermite function of index v. 

We have related the Hermite functions to the functions 



I(v, a,X)= x 



„v—\—ax—\x* 



dx = (2\)-" /2 T{v)h- v 



in Carr, Gcman, Madan and Yor (2005). 
We may therefore write 



h- Y (V2\) = (2A)* ^/(Y^A.A) 
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It follows that 

dx + x) 



-i—A -XX 
2 2 g AA 



It follows that 



, xi (A) = 2 r A- r ^ + 5) 



r(y)r(|) 



ffr(y,2A,A) 



We therefore evaluate 



E 



B 2 '' — 



r(F)r(i; 



I[Y,B'y, 



B 2 y 



(13) 



Putting together the result of equation (|13|l and equation (|12fl we get the 
results for the CGMY subordinator J5J. 



4 Simulating CGMY using Rosinski Rejection 

We suppose that we have two Levy measures Q(dx), Qo(dx) with the property 
that 

dO 

dQa 

and this is our case, then it is shown in Rosinski that we may simulate the paths 
of Q from those of Qq by only accepting all jumps x in the paths of Qq for which 

dQ 

(X) > w 



dQo 

where w is an independent draw from a uniform distribution. 
For our case we have that 

dV\ „ r „7i 

— - = E \e yZ ] < 1 
dv 

and so accept all jumps in the paths of vq for which 

E [e- yZ ] > w 

The detailed algorithm is for parameters C, G, M, Y to first define the time 
step to be C, 

t = C. 
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Then we let 

A = 

= G + M_ 
2 

We next simulate at time t from the one-sided stable subordinator with Levy 
measure 

dy 



V 2 + 

For this we let e — .0001 and truncate jumps below e replacing them by 
their expected value at a rate of 



r 1 
d = / y^zr d y 

J0 !)! + 



1- ^ 

1 2 



For the arrival rate of jumps we have an arrival rate A of 

,oo 1 

x = / -^T7 dy 

Je y 2 

2 1 

Y e\ 

The interval jump times are exponential and are simulated by 

U = -^log(l - U 2 i) 

for an independent uniform sequence U2i- The actual jump times are 

i 



For the jump magnitude we simulate from the normalized Levy measure the 
jump size yj given by 

e 

Vi = ~ rx 

for an independent uniform sequence «y. 

The process S(t) for the stable subordinator is given by 

oo 

S{t) = dt + J2vj 1 r j <t 
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We now get the CGMY subordinator H(t) by 

OO 

H(t) = dt + ^ y:j l T]<t l h(y)>U3o 

w r(y)r(i) v 2 J v 2 

for an independent uniform sequence u 3 j 

Finally we simulate the CGMY random variable by 

X = AH(t) + ^H{t)z 

for a draw z of a standard normal random variable. 



5 The Meixner Process as a Time Changed Brow- 
nian Motion 

We consider the Meixner Process (Schoutens and Teugels (1998), Pitman and 
Yor (2003)) as a time changed Brownian motion. The Levy measure of the 
Meixner process is 

, „ exp (-x) 
k{x) =5- ' ••' ' 



xsinh(^) 
The characteristic function is given by 

4>Mei X neM = E[e^} 



cos(&/2) 



cosh(au — ib)/2) 



To see this process as a time changed Brownian motion we wish to identify 
l(u) the Levy measure of a subordinator such that 

M» = / -^=<>sv[- {x ~ Ay) )l(y)dy 



r 1 ( 

I —= exp - 

J-oo V^V \ 



2y 



^ 1 / ^A-i-^ ] ' [ ' lUh < 



Hence we set 

A = - 

a 



and seek to write 

1 f°° 1 fx 2 A 2 y\ , s 
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We transform the left hand side of 114|) as follows. 
We recall that 



Cx 



where T. 



(3) 



c 



inf ji| 



Then we write 



sinh(Cir) 



= E 



exp 



2 c 



for i?^ 3) 



the BES(Z) process. 



1 

:sinh(2^) 



(^) sinh (*2 
(5a 1 



7T Or 

5a 1 

7T X 2 



exp 



exp 



with C = — . Denote by 9{h)dh the law of . We may then write 



" 2 Jc 



; sinh I 



where 



We recall that 



5a 



du 



— I "7T ex P 



exp 



x 2 C 



(3) 



n r°° r°° / r 2 

-/ / 6>(*)diexp - — (u + C 2 t) 
K Jo Jo V 2 



— / du£ 

2tt 

— / du 

2' 

(5(7 

2^ 

(5(7 

2^ 

(5(7 

2^ 



exp(-y (u + CTJ 



2/p(3) 



.-/'/ / ^exp 



duexp ( - — 



2 

du 
C 2 



C 2 

v — u 



C' 2 



duexp ( - — 



dh9(h) 



duexp ( — ) 



5a 
27 

-P ( T, 

2tt 

5a 
2^' 



6(h)dh 

(3) < M 



a 2 / 



>(3) 



j,(3) (to) 



( max t <i i? 



(3) 
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We now transform the right hand side of (|14|) to write 
f°° 1 fx 2 A 2 y\ J ,^ 1 r 1 / x 2 v A 2 \/l 

From the uniqueness of Laplace transforms we deduce that 
?( w ) = _Uexp(^W- 



dv 



or 



l W = \ -^V - exp 



u 

2ir 8a „ / , ,c 3 i ^ r-\ ( A 2 u 
--p(Mf)>C^)exp(^- — 



exp 



where 



V27TU 3 

r (3) - ^ / A 2 u 



g(u) = P Ml 6 ' > CVu) exp 



For the absolute continuity of our subordinator with respect to the one sided 
stable i subordinator we require that 



I , \ 2 

f= [Vgiu) - lj du < 



For this we observe that 

2 



R-l) < |5(«)-1| 
= l-ff(u) 

= l-p(Mf)>CVS)exp / S "" 



2 

= 1 - exp ( ) + exp [-^) (i-P ( M i aJ > CV^ 

The first part is clearly integrable with respect to (^72-)and for the second we 
observe that as 

X k P (T > A) <E [T k ] 

that 

P\ ^ > -4- ) = P (if 3 > -p^—] < Ku\ for all k 
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For the simulation of Meixner as a time changed Brownian motion we would 
wish to evaluate 

p (jifW > CV5) = p(^^<7^-) 

\ ) V( M i ) 2 Cu ) 

- p (^^m 

= 1^£> 

CO 

— oo 

For the last equality we refer to Pitman and Yor (2003). 

6 Simulation of the Meixner Process 

The simulation strategy is similar to that employed in section 3 for CGMY, 
except that here we simulate first the jumps of the one sided stable | with Levy 
density 

, . . 5a 

k{x) = , x > 0. 

V 27TX 3 

We approximate the small jumps of the subordinator using the drift 

The arrival rate for the jumps above e is 

A = Sa\ — 



V 7T£ 

and the jump sizes for the one sided stable(|) are 

for an independent uniform sequence Uj. 

We then evaluate the function g{y) at the point yj and define the time change 
variable 

r = S + ^2vjlg(vi)>w j 

3 

for yet another independent uniform sequence Wj. We note that the function 
g(y) only use the parameters a, b and is independent of the parameter d. 
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The value of the Meixner random variable or equivalently the unit time level 
of the process is then generated as 

b r- 

X = -t + ^Fz 
a 

where z is an independent standard normal variate. 

7 Results of Simulations 

For both the CGMY and Meixner processes we present in this section the 
results of simulating the processes at typical parameter values obtained on 
calibrating option prices on the S&P 500 index. The parameter values for the 
CGMY are C = 1, G = 5, M = 10, and Y = .5. The parameters for the 
Meixner were a = .25, b = —1.5 and 5 = 1. 

We present graphs ill 12^ for a weekly time step h = .02 of the simulated and 
actual densities as well as chi square tests of the hypothesis that the sample was 
drawn from the respective densities. The solid lines are the theoretical density 
while the data points are indicated by dots. The sample sizes in both cases were 
5000. The range for both the CGMY and Meixner returns was 25%. In both 
cases we used 100 cells and employed those with more than five observations for 
the test. The CGMY had a chisquare statistic of 42.0122 with 56 degrees of 
freedom and a p — value of .9172. For the Meixner the test statistic was 78.70 
with 84 degrees of freedom and a p — value of .6427. 
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Figure 1: CGMY simulation as time changed Brownian Motion using shaved 
one sided stable Y/2. 
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Figure 2: Meixner simulation as time changed Brownian motion using shaved 
one sided stable 1/2. 
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